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ABSTRACT 

The very small braking index of PSR J1734-3333, n = 0.9 ± 0.2, challenges 
the current braking mechanisms of pulsars. We present a possible interpretation 
that this pulsar is surrounded by a fall-back disk and braked by the disk. We 
analyzed the interaction between the neutron star magnetosphere and a possible 
fall-back disk, and found that the braking index of the neutron star is very small 
in the middle stage of propeller phase. In this regime, the braking index of PSR 
J 1734-3333 will increase in the later stage of propeller phase, when it evolves to 
the region of anomalous X-ray pulsars and soft gamma repeaters in the P — P 
diagram. The mass of the disk around PSR J1734-3333 in our model is about 
20 - 30M e , similar to the observed disk mass around AXP 4U 0142+61. 

Subject headings: pulsars: individual (PSR J1734-3333) - stars: evolution - stars: 
magnetars - stars: neutron 
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Introduction 



PSR J1734-3333 is a rotation-powered radio pulsar with a period P = 1.17 s and 
period derivative P = 2.28 x 10~ 12 . It is located between the normal rotation-powered 
radio pulsars and anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs) in 
the P — P diagram. PSR J1734-3333 has been observed regularly since 1997 by using 
the 64 m telescope at Parkes and the 76 m telescope at Jodrell Bank. No glitch occurs 
during these years, which ensures the accuracy of measuring z> and u, where v = 1/P is the 
rotation frequency of the pulsar. The braking index of a pulsar, n = vvjv 2 ^ describes the 
dependence of the braking torque on r otation frequency. Fo r PSR J1734-3333, n = 0.9 ±0.2, 



based on 13.5 years observation data (jEspinoza et al. 



201lh . 



Because of the observational 



braking indices (Lvne et al. 



Weltevrede et al. 



1993 



difficu 



1996 



ties, only seven pulsars have reliable pub 



Middleditch et al. 



2006 



Livingstone et al. 



ished 



2007 



20111 ) before PSR J1734-3333. All of th em have a value between 1 and 



3 except PSR J0537-6910, which has glitched frequently ( IMiddleditch et al. 



20061 ). The 



braking index value of PSR J1734-3333 is too small compared to that from spin-down 
models and the measured values of other pulsars. The pure magnetic braking with a 



Xu fc Qiaol ( 1200 ll ) considered the 



constant dipolar magnetic field should have n = 3. 
braking torques due to magnetodipole radiation and the unipolar generator, and predicted 
1 < n < 3. Other factors may also affect the braking index, such as the stellar winds, 



systematic variations o: 



( [Manchester et al. 



1985 



;he moment of inertia or m agnetic field, and non-dipolar braking 



Blandford fc Romani 



19881 ). Therefore, the small braking index of 



PSR J1734-3333 is a challenge to the current braking mechanisms of pulsars. 

From the definition of the braking index, P = kP 2 ~ n . When a pulsar evolves, its 
location in the P — P diagram will move according to n. The slope of the evolutionary 
path in the diagram is 2 — If PSR J1734-3333 keeps the present braking index value, it 



-4- 



woulc 



fl201lh sug 



1995 



dramatically evolve into the AXPs & SGRs region in ~ 20, 00 yr. 



gested that PSR J1734-3333 may be a potential magnetar ([Thompson fc Duncan 



Espinoza et al 



1996), whose magnetic field has been buried under the surface due to large accretion 
shortly after the supernova explosion, and is relaxing out of the surface at present. The 
increasing of magnetic field strength could also result in the small braking index value of 
0.9. In this paper, we propose another interpretation for the small braking index that there 
may be a fall-back disk around PSR J1734-3333 which brakes the pulsar. 



2. Braking PSR J 1734-3333 by a possible fall-back disk 

2.1. The fall-back disk 

It is reasonable to assume that the mass could not be perfectly ejected during the 
supernova explosion and that a small amount of mass can fall back. A part of the fall-back 
material carry a sufficient angular momentum and can rotate around the young neutron 



ill-back disk is unclear, bu 


t could be larger 


( Chevalier 


1989; 


Lin et al. 


1991; 


Wang et al. 



20061 ). PSR J1734-3333 is a young pulsar, with a characteristic age r c ~ 8000 yr. There 



may be a fall-back disk around it. 

The fall-back disks are very difficult to be detected due to their extremely weak 
emission. About half of the AXPs SGRs have infrared /optical coun t erparts obse r ved 



(I Wang fc Chakrabarty 



2002 



Israel et al 



2004 



2005 



Mereghetti 



2011a). 



Wang et al 



(120061 ) 



observed the mid-infrared emission from a debris disk around AXP 4U 0142+61, the 



brightest known AXP. The disk mass is of the order of ~ 10M®. iKaplan et al. 



J2009|) 



detected the counterpart to the AXP IE 2259+586 at 4.5 jum, which can be explained by 
the passive X-ray- heated dust disk model developed for 4U 0142+61. On the theoretical 
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point of view, the existence of disks around AXPs & SGRs can help to understand their 
distribution on the P — P diagram, e.g. the origin of SGR 0418+5729 with a dipolar 



magn etic field less than 7.5 x 10 G ( Ivan der Horst et al 



2010; 



Rea et al. 



20101 : 



Alpar et al. 



201ll ). It is possible that PRS J1734-3333 will evolve to be a AXP or SGR, and now has a 
not-yet-detected disk around it. 



When a fall-back disk i s formed around a pulsar, the accretion rate will evolve with 



time. 



Cannizzo et al. 



( 119901 ) found that the accretion rate declines self-similarly according 
to M oc t~ a , where a > 1, when the disk is under the influence of v i scous processes. We 
parameterize the accretion rate in the same way as 



Chatterj ee et al 



((2000), i.e. 



M = M , 
M = M (t/T)- 



< t < T, 



(1) 



a , t > T, 



where T is the dynamical time in the inner parts of the disk, and M is a constant accretion 
rate. With these assumptions, the total disk mass is 



and the residual to accrete mass is 



Mdt 



-M T, 



(2) 



M d (t) 



Mdt 

ft — 1 



(3) 



In this model, the super-Eddington accretion could occur at early accretion phase. The 
initial parameters Mo and T are coupling, and they are insensitive to the results at a late 



stage. 



Cannizzo et al. 



( 119901 ) found a = 19/16 and 1.25 for a disk in which the opacity 
is dominated by electron scattering and for a Kramers opacity, respectively. We choose 
a = 7/6, which i s very nearly the above values and could simply the calculations, as 



Chatterj ee et al. 



(120001 ) has shown. Our calculations show that small changes of the a value 



do not effect the results significantly. 
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2.2. Period evolution of PSR J1734-3333 with a fall-back disk 

The fall-back disk will interact with the pulsar magnetosphere. The influence on the 
spin evolution of a pulsar is 

m = j, (4) 

where I and Cl are the moment of inertia and angular velocity derivative of the pulsar, and 
J is the torque acting on the pulsar by the disk. To study the spin evolution of a pulsar in 
detail, we introduce two important radii: one is the corotate radius R co where the corotate 
velocity is equal to the Kepler velocity, and the other is the magnetosphere radius R m inside 
which the accretion disk cannot exist. The corotate radius 

R co = (GM/4tt 2 ) 1/3 P 2/3 « 1.7 x 10 8 P 2/3 cm, (5) 

where G is the gravitational constant, M and P are the mass and period of the star. The 
magnetosphere radius 

R m « 0.5r A « 8 x 10 7 )i% 7 M~ 2/? cm (6) 



( iFrank et al.lll992l ). where r& is the Alfven radius inside which the flow is dominated by the 
magnetic field, /130 the magnetic moment of the star in units of 10 30 G cm 3 (for a dipole 
magnetic field, fi = BR 3 ), and Mig the accretion rate in units of 10 18 g/s. 

The sizes of R co and R m will change with the evolution of the pulsar and the disk, 
thereby the accretion process would undergo three phases successively: the accretion phase, 
the propeller phase and the tracking phase. At the early stage, R m < R co , the accretion 
flow can fall onto the star surface, and the gravitational energy is released. At the initial 
stage of the accretion phase, M is larger than the Eddington limit ( Me ~ 10 18 g/s), so that 
only part of the accretion flow Me can reach the star surface, and the excessive material 
is blown away by the radiation pressure from the star. The accreted material transfers its 
angular momentum to the neutron star, while the blown away material may take away 
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angular momentum of the neutron star. This process has been discussed by iLiu et al. 



( 120121 ). We found that the period changes very little at this stage because the time scale for 
such process is very short. The magnetosphere radius R m increases faster than the corotate 
radius R co , and would exceed R co after a short accretion phase. When R m > R co , the whole 
accretion flow c annot reach the star surface , and the material would be ejected out by the 
propeller effect (jlllarionov fc Sunyaevlll975l ). The propeller effect can brake the star to a 
relatively long spin period (typical several seconds), until R m approximates to R co and the 
propeller torque fades away. i? m could never equal to R co , but changes with R co . This is 
the tracking phase. In the whole evolution process, the propeller phase is the most interest 
and important stage to brake the star. 

The propeller effect betwe en a rotating magnet i c neu tron star and the accreting 



material was first proposed by 



Illarionov fc Sunyaevl (119751 ). to explain why the number 



of galactic X-ray sources is much less than the expected number of neutron stars and 
black holes in binary systems. They suggested that most of neutron star binary systems 
are in propeller phase, so that the accreting material is thrown away from the neutron 
star by t he rotating magnetosphe re. They could not emit enough X-ray luminosity to be 



detected. 
B1259-63 in a 
is yet unclear. 



Manchester et al. 



( 119951 ) found evidence for a propeller-torque spin-down of PSR 
lighly eccent r ic orb it bi nary system. However, t he value of propeller-torque 



Menou et al 



(1199a ) and 



Chatterj ee et al 



( 120001 ) proposed a simple idealized 



formula for the torque, that 



j = 2MR 2 m Q K (R m )[l 



n 



(7) 



where ^(flm) is the Kepler angular velocity at R m , and Q the angular velocity of the 
star. In this formula the torque is equivalent to that for an ideal process. The accretion 
material rotates with f2x before it drops to R m , and corotates with the star when it drops 
to i? m , and finally is thrown away from the neutron star. The torque is zero when the star 
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angular velocity Q equals to ^(-Rm), and would increase with the difference between Q 
and Q}<i{Rm)- This is reasonable and feasible. However, this torque seems to be too large 
to be true when Q » Q^(R m ), considering the facts that the magnetic energy density is 
comparable to the kinetic energy density of free-fall material at R m , which is two times the 
Kepler motion energy density. According to this formula, the kinetic energy density would 
become much larger than the magnetic energy density when Q » Q^(R m ), which seems to 
be more or less unnatural. Considering the reasonable physics and unusually large torque 
for Q » fi K (-Rm), we modify the formula for the propeller torque as 

j = 2MR 2 Jl K {R m ){\ - [o-^r]*} (8) 

with an additional parameter x, < \ < 1. The torque can be adjusted by the parameter 
X- 

Given the accretion history from equation (1) and the torque formula (8), the spin 
evolution of a pulsar can be calculated using equations (4) and (6). The magnetic dipole 
radiation torque is insignificant compared with the propeller torque, especially when it 
spins slowly. The initial accretion phase is very short (~1 year), and its effect on the period 
evolution is very small. We ignore the magnetic dipole radiation and the initial accretion 
phase in the following calculations. 

Figure 1 shows the evolution paths of a neutron star surrounded by a fall-back disk in 
the P — P diagram, with initial parameters Pq = 15 ms, T = 1 s, and a = 7/6. We tried 
X = 1/3 and x — 1/2, whose paths can go across the location o f PSR J1734-3333, an d the 



braking indexes coincide with t 



re observed value n = 0.9 ± 0.2 (jEspinoza et a 



comp arison, the case of x — 1 (jChatterjee et al 



2000|) and the case of x — 2 (lErtan et al 



2011"). For 



20091 ) are also calculated and shown in Figure 1, in which the braking indexes can not be 
equal to the observed value of PSR J1734-3333. As Figure 1 shows, the propeller phase 
can be divided into three stages according to the slope of the evolution path in the P — P 
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diagram. In the first stage, the path goes down quickly until the slope increases from very 
negative to zero. In the second stage the path goes up, where the slope increases from zero 
to a value larger than 1, i.e. n < 1, and decreases to zero at the end. In the last stage the 



path goes down 



again, and the s 



ope k eeps decreasing, until the neutron star evolves to the 



20001 ) . The fitting parameters and main results are given 



tracking phase (jChatterjee et al. 
in Table 1. 

Figure 2 shows the period evolution with the different parameters given in Table 1. 
The period of a neutron star surrounded by a fall-back disk could evolve to several seconds 
in about tens of thousands of years and could stay in the period range for a very long time. 



3. Discussions and Conclusions 



In this paper we proposed a moderate formula for the propeller torque considering 
the energy dens ities of the magnetosphere and the accre t ion materia l , motivated 



previous works (IMenou et al. 



1999 



Chatterj ee et al. 



2000 



Alpar 



2001 



, motivated 


by 


Ertan et al. 


2009) 



A suggested fall-back disk around PSR J1734-3333 with such a torque can explain the 
observed small braking index. In most of the second stage of the propeller phase, the 
braking index is about 1 as shown in Figure 1. In the third stage of the propeller phase, 
the neutron star can evolve into the region of AXPs & SGRs in the P — P diagram, and 
its period derivative can change in a large range. This evolution scenario can explain the 
periods clustering of AXPs & SGRs in the P — P diagram. 

As shown in Ta ble 1, the initial disk masses are smaller than O.1M which agrees with 



Chevalierl ( 119891 ) and 



Lin et al 



f ll99lh . When PSR J1734-3333 evolves into the AXP & SGR 



region in about 20,000 years later, the disk mass would be of t 



re order of 10M ffi , similar to 



the observed disk mass around AXP 4U 0142+61 (jWang et al. 



20061 ). The propeller torque 
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of a fall-back disk can modify the period derivative, which makes the apparent dipole 
magnetic field strength much stronger than the real field strength. Our model implies that 
the magnetic field strength of PSR J1734-3333 is only about 10 12 G, much smaller than the 
expected value from the pure magnetic braking. 
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supported by the National Basic Research Program of China (Grant Nos. 2009CB824800, 
2012CB821800), the National Natural Science Foundation of China (Grant Nos. 11225314, 
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Fig. 1. — Evolution paths of a neutron star surrounded by a fall-back disk with different 
propeller torques, given initial parameters Pq = 15 ms, T = 1 s, and a = 7/6. To cross 
the location of PSR Jl 734-3333 in the diagram, different \ needs different disk mass and 
magnetic field strength. The fitting parameters and main results are given in Table 1. The 
two short lines near the star indicate the observational upper and lower limits of the braking 
index of PS R J1734-3333. The i nform ation of radio pulsars was taken from the ATNF pulsar 



catalogue (Manchester et al 



( 120051 ). http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ ), 



data of AXPs and SGRs from the McGill SGR/AXP Online Catalog 
(http: / /www.physic s.mcKill.ca/ p ulsar/m ag;netar/main.html), and X-ray Dim Isolated 



Neutron Stars from 



Mereghettil (1201 lbf ) 



-14- 



O 10 

CD 




10 10 
time (year) 



Fig. 2. — Period evolution of a neutron star with P = 15 ms surrounded by a disk with 
different parameters given in Table 1. 



Table 1: Fitting parameters and main results. 



X 


1^30 


M d (M ) 


t(kyr) 


n(t) 


M d (t)(M e ) 


1/3 


3.3 


0.035 


59.8 


0.92 


28 


1/2 


1.0 


0.025 


50.8 


1.07 


20 


1 


0.48 


0.0074 


29.0 


1.57 


8 


2 


0.36 


0.014 


10.7 


2.8 


15 



Note. — Md, t, n(t) and M^(t) are the initial mass of fall-back disk, age, current braking index, and 
residual disk mass, respectively. 



